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ASYMPTOTIC DIRECTIONS, MONGE-AMPERE EQUATIONS 
AND THE GEOMETRY OF DIFFEOMORPHISM GROUPS 

BORIS KHESIN AND GERARD MISIOLEK 

Abstract. In this note we obtain the characterization for asymptotic directions 
on various subgroups of the diffeomorphism group. We give a simple proof of non- 
existence of such directions for area-preserving diffeomorphisms of closed surfaces 
of non-zero curvature. Finally, we exhibit the common origin of the Monge- Ampere 
equations in 2D fluid dynamics and mass transport. 



I. Asymptotic directions on diffeomorphism groups 

The group of volume-preserving diffeomorphisms of a Riemannian manifold plays 
a fundamental role in the geometrical study of the Euler equation of hydrodynamics 
on the manifold |Arj . In this paper we consider another equation, the Monge- Ampere 
equation, and discuss its universality in the context of diffeomorphism groups. This 
equation occurs in two main contexts: as the equation of asymptotic directions in 
2D hydrodynamics and in the optimal transport problem in any dimension. 

1.1. Characterization of asymptotic directions on subgroups. Let M be a 

compact n-dimensional manifold without boundary and equipped with a Riemann- 
ian metric g. Let T>(M) be (the connected component of the identity in) the group 
of all diffeomorphisms of M. Its tangent space at the identity diffeomorphism con- 
sists of smooth vector fields on M. The tangent space at a point rj consists of vector 
fields "reparameterized by 77," i.e. of the maps X v : M — ► TM with X n (x) G T V ^M. 
Define the "flat" L 2 -metric on V(M) by assigning to each tangent space the inner 
product 



1.1) g v (X v ,Y v ):= [ g v{x) (X v (x),Y v (x)) dV(x) 

J M 



where dV denotes the Riemannian volume form on M. 

Let ST>(M) be the subgroup of volume-preserving diffeomorphisms. The restric- 
tion of the L 2 -metric to this subgroup is right-invariant and of particular importance 
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in hydrodynamics. In jArj Arnold showed that geodesies in ST>(M) correspond to 
motions of an ideal fluid in M described by the Euler equations 

d t X + V X X = - Vp, div X = 

for the fluid velocity field X. 

As shown by Ebin and Marsden in JEM] the group T> S (M) of all diffeomorphisms 
of Sobolev class H s (as well as its various subgroups) can be viewed for s > n/2 + 1 as 
an infinite dimensional Hilbert manifold. All of the arguments in this paper can be 
rigorously developed in the Sobolev framework. However, to present the geometric 
ideas we will keep things formal and drop the index s in what follows. 

While the Euler equation depends only on the intrinsic Riemannian geometry of 
ST>(M), it is also of interest to study its exterior geometry as a Riemannian sub- 
manifold in T>(M). In particular, one can consider asymptotic directions in ST>(M). 
A vector tangent to a Riemannian submanifold is asymptotic if the geodesies issued 
in the direction of this vector, one in the submanifold and the other in the ambi- 
ent manifold, have a second order of tangency. (Note that in general two geodesies 
with a common tangent will have only a simple, i.e. first order, tangency.) More 
formally, asymptotic vectors are singled out by the condition that the second funda- 
mental form evaluated on these vectors is zero. A curve whose tangent is asymptotic 
at each point is called an asymptotic line. An asymptotic line is a geodesic in the 
submanifold if and only if it is also a geodesic in the ambient manifold. 

A description of asymptotic vectors in the group of volume-preserving diffeomor- 
phisms was given by Bao and Ratiu in |BRj . while asymptotic geodesies in ST>(M) 
were studied in |Mlj (as pressure-constant flows). We begin with the following 
convenient characterization of these vectors. 

Theorem 1. [BRJ A vector field X on a manifold M is an asymptotic direction for 
ST>(M) at the identity diffeomorphism if and only if 

(1.2) divVxA" = divX = 0. 

Similarly, a vector X n is asymptotic to ST>(M) at the diffeomorphism r\ if and 
only if its right translation X := X v o rj' 1 is asymptotic to ST>(M) at the identity. 
If M is two-dimensional then more can be said: 

Corollary 2. [IBRj The stream function ip of an asymptotic vector field X in 2D 
satisfies a Monge- Ampere equation 

(1.3) det[DV] = ^|W| 2 , 

where g = det(gij), K is the Gaussian curvature function of M and det[D 2 ip] is 
the Hessian of ip for the field X = JVip with respect to the symplectic area form 
oj = dV on M. 

Below we give a characterization of asymptotic vectors in a general setting. Let 
B be a subgroup of V(M) and let b denote its Lie algebra. We assume that b is a 
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closed subspace, and therefore it has an orthogonal complement with respect to the 
L 2 inner product (jl.lj) . 

Theorem 3. A vector field X on M is asymptotic for B if and only if 

(1.4) Xeb and V x Xeb. 

If B = ST>(M) then the Lie algebra b consists of divergence-free vector fields and 
we recover Theorem ^ Moreover, if M has a boundary then the diffeomorphisms 
from the subgroup B leave it invariant and hence the fields from the subalgebra b 
are tangent to the boundary of M. 

Corollary 4. |BRj If M has a boundary dM , then X is asymptotic for ST>(M) if 
in addition to equations Al.fy) X satisfies the conditions 

(1.5) g(X,n)=g(V x X,n) = 
where n is the normal to dM. 

Suppose next that M is a symplectic manifold of dimension 2n equipped with a 
symplectic 2-form u. 

Corollary 5. a) Let B = Symp(M) be the subgroup of symplectic diffeomorphisms 
with the corresponding Lie subalgebra b = symp(M) of vector fields preserving u. 
Then a vector field X is asymptotic if and only if 

Lxoj = and L^ x xoj = 0. 

b) Let B = Ham(M) be the (generally speaking, smaller) subgroup of Hamiltonian 
diffeomorphisms. Then a vector field X on M is asymptotic if and only if 

uj(X, •) = dip and u(VxX, •) = dip 

for some smooth functions ip and <fi on M . 

If dim M = 2 any divergence-free field is locally Hamiltonian and the latter equa- 
tion is rewritten in (jl.3|) as the Monge-Ampere equation on its Hamiltonian, or 
stream function. 

Curiously, this Monge-Ampere property does not survive when passing to Hamil- 
tonian fields in higher dimensions, unlike what was conjectured at the end of |AKj . 
Indeed, already in the flat 4-dimensional case one obtains a system of three equations 
on pairwise products of second derivatives, while the corresponding Monge-Ampere 
with det [-D 2 ^] would include the four-term products of second derivatives. 

Another interesting example is provided by a contact manifold M with a contact 
structure (i.e. maximally non-integrable distribution of hyperplanes) r. 

Corollary 6. Let (M, r) be a contact manifold, and B a subgroup consisting of 
contact diffeomorphisms. Then the vector field X is asymptotic to B if and only if 

both X and V xX are contact. 
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If a denotes a 1-form defining the contact structure r then the Lie algebra b := 
{X : Lx<y = fot}. For a contact vector field X consider its contact Hamiltonian 
function Kx '■= a(X) with respect to the form a. By rewriting the condition 
VyX e b for asymptotic vectors on B one can obtain an analog of the Monge- 
Ampere equation (jl.3j) on the Hamiltonian in the contact 3-dimensional case. 

Proof of Theorem OH The Lie algebra of T>(M) splits orthogonally into b © b -1 
inducing a corresponding splitting at any point in B by right translations. Let V 
denote the smooth Levi-Civita connection of the L 2 -metric (jl.lj) on T>(M) (see Ebin 
and Marsden [EMJ). Let X and Y be two elements of b and extend them to right- 
invariant vector fields X v = X o rj and Y v — Y o 77 on B. Decomposing into unique 
tangential and normal components (Gauss equation) we obtain 

V Xv Y v = P v (V Xv Y v ) + (Vx/,) 1 . 

where P v is the projection onto the tangent space to B at the point 77. The tangential 
component can be used to define a smooth right-invariant connection on B while 
the normal component (Vx l^) 1- is interpreted as the second fundamental form of 
B in V(M). Since 

V x Y v = (V x Y)or ] 

we see that at the identity the L 2 covariant derivative V is given by VxY. Recall 
that a vector X in b is asymptotic if the second fundamental form evaluated at X is 
zero. This implies that (VxX) 1 - = and Theorem El follows. □ 

Remark 7. One can see how the above general consideration works in the case 
B = ST>(M). First, recall that an arbitrary vector field X on M can be decomposed 
into L 2 -orthogonal divergence-free and gradient parts 

X = P ld (X)+X v 

where Pid is now the projection onto the divergence-free fields, and X v is the gra- 
dient part of X: explicitly X v = VA _1 divX. By right invariance this induces 
a corresponding splitting at any point in ST>(M). The tangential component of 
V defines a smooth right-invariant connection on ST>(M). Its normal component 
is the second fundamental form of the subgroup ST>(M) in T>(M). Thus, for any 
right-invariant vector fields X v = X o rj and Y v = Y o 77 on ST>(M) we again have 
the Gauss equation 

(1-6) V X Y V = P V (V X Y V ) + s(X v , Y v ) 

where 

V Xr ,Y ri =(VxY)o V and s(X v , Y„) = (V X Y) V o 77. 

Since a vector X tangent at the identity r] — id is asymptotic if s(X,X) = we 
immediately obtain the characterization in Theorem ^ 
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Remark 8. Asymptotic geodesies (or pressure-constant flows) are of interest in the 
Lagrangian approach to hydrodynamic stability theory. Typical examples are plane- 
parallel flows on the flat torus T 2 . These flows can be considered unstable in the 
following sense. For any such flow it can be shown that sectional curvatures of 
SV(T 2 ) along the corresponding geodesic are always non-positive. Therefore, by 
a suitable variant of the Rauch comparison argument, all linear perturbations in 
Lagrangian coordinates (Jacobi fields along the geodesic) must grow at least linearly 
in time. On the other hand, there are flows with non-positive curvatures for which 
one can show that the growth must be at most polynomial (see Preston |Prj for 
details). 

1.2. Non-existence of asymptotic directions. Asymptotic directions are not 
always in good supply. In BRJ it is shown that if M is a two-dimensional compact 
surface of revolution without boundary then any axially symmetric smooth solution 
of the Monge- Ampere equation is constant away from the cylindrical (i.e. fixed ra- 
dius) bands of the surface. The strongest result in this direction so far states that for 
a compact closed surface M of positive curvature there are no asymptotic directions 
on ST>(M), see jPaj . A similar result holds for surfaces with boundary. (Palmer also 
showed that there is no direct analog of this result in higher dimensions: e.g. every 
left-invariant vector field on a compact Lie group equipped with a bi-invariant met- 
ric solves the equations (|1.2|) . i.e. is asymptotic. This had been observed previously 
for the three sphere S 3 in fHl]-) 

Here we prove the following generalization of Palmer's result. 

Theorem 9. If M is a compact closed surface of nowhere zero curvature K , then 
the Monge-Ampere equation admits no non-constant solutions. In particular, 

in this case ST>(M) has no asymptotic directions. 

Note that the surface can be of any non-zero Euler characteristic (the case of the 
torus is ruled out by the Gauss-Bonnet theorem). 

Proof. Recall that for any vector field X on M we have the identity 

(1.7) dwV x X = r(X,X) +tr(DX) 2 + L x (dwX) 

where r denotes the Ricci curvature of the metric g and Lx is the Lie derivative 
along X (see for example Taj). If X is divergence-free then the last term on the 
right side of ()1.7|) drops out. 

Consider the length function / := g(X,X). Since M is compact / must attain 
a maximum at some point Xq. Choosing normal coordinates at that point with 
9ij( x o) — $ij an d dgij/dx k (x ) = 0, we obtain 

= df(x ) = 2 £ X\x ) ^r(xo) dx k . 

jk 

Since at the point where / has a maximum we must also have X(xq) ^ 0, this 
implies that the Jacobi matrix DX is degenerate at xq. Therefore, rearranging the 
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terms of the trace and using the fact that the divergence of X is zero, we get 

r) Xi r) X* 

tr(DX) 2 (x ) = J2 = - 2det [ M (^o)] = 0. 

ij 

Substituting into the relation ()1.7|) we find that 

div VxX(xq) = K(x ) g(X, X) (x ) , 

because in two dimensions the Ricci and Gaussian curvatures coincide. However, if 
X is asymptotic this implies that 

= K(x o )g(X J X)(x o ) 

contradicting the assumption K ^ on M, and in particular at the point Xq. □ 

Remark 10. Note that the proof above works for a C 1 field X, thus improving on 
the C 2 assumption used in jPaj. 

The following is a similar result for the case of a manifold with boundary: 

Theorem 11. Let M be a compact surface of nowhere zero curvature K with smooth 
boundary. Assume that the geodesic curvature k g of dM vanishes, at most, at finitely 
many points. Then the Monge- Ampere equation M.^l along with the boundary con- 
dition / li.,5]) admits no non-constant solutions. 

Proof. The beginning of the proof follows |Paj . Assuming that such an X exists we 
show that X has to vanish on the boundary. Indeed, for any point x G dM where 
X(x) ^ we have that X is tangent to dM and 

= (V x X,n)(x) = k g (x) ■ g(X,X)(x) . 

The latter shows that the geodesic curvature k g must vanish at x. By the assumption, 
this implies that X can be non-zero only at finitely many points of the boundary. 
Hence it is identically zero on dM by continuity. 

The proof of the boundary-free case is now applicable: the function / := g(X, X) 
is also zero on the boundary dM and must therefore attain a maximum in the 
interior of M. □ 

2. Mass transport and foliations on diffeomorphism groups 

2.1. Monge-Ampere equation in optimal mass transport. A somewhat dif- 
ferent Monge- Ampere equation arises in the theory of mass transport. Namely, let 
dfi(x) = m(x)dV(x) and du(y) = n(y)dV(y) be two smooth measures on a manifold 
M of any dimension. The Jacobian of a map rj which sends the measure /i to the 
other v = f]*[i satisfies the relation 

(2.1) n(rj(x)) det[Dr](x)] = m{x). 

The property of optimality means that the 1-parameter family of maps t] t (x) de- 
scribes a geodesic curve on the space of densities with respect to the Wasserstein 
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L 2 -metric. The latter is the 'transport' L 2 metric on the space of densities: the 
distance between two densities on M is the cost of transporting one of them to the 
other with the L 2 -cost function. 

If M is a domain in M. k one can see that an optimal map 77 : M — > M has to be 
the gradient rj = V0 of a convex function <p, i.e. the equation on the potential 4> 
assumes the Monge- Ampere form (see, e.g, [BrJ): 



(2.2) det[D 2 0] = ZJ^ - 
v 1 11 n(V0(x)) 

In the case of an arbitrary manifold M the optimality condition is expressed in 
terms of convexity of the function with respect to the metric on M, see |Mcj . The 
corresponding Monge- Ampere equation has the same form as above. 

2.2. Universality of the (pre-) Monge- Ampere equation. Although the two 
Monge-Ampere equations discussed above look rather different (equation (J1.3j) is 
for asymptotic vector fields, while equation (J2.2)) is for potentials of optimal diffeo- 
morphisms), both have a common origin. It turns out that they can be viewed as 
projections of the dispersionless Burgers equation: 

(2.3) d t X + V X X = , 

which can be thought of as a pre-Monge-Ampere equation. This equation describes 
geodesies with respect to the "flat" L 2 -metric (II. lj) on the group of all diffeomor- 
phisms T>(M) for M of any dimension. According to this equation each fluid particle 
moves along a geodesic in M. 

Note that there is a natural fibration on T>(M) given by the projection n onto 
densities V(M). Two diffeomorphisms belong to the same fiber if they move a fixed 
density (say, the constant density O for a compact M) to one and the same density. 
In particular, for a compact M one considers densities with the same total mass, and 
<SX>(M) = 7r _1 (0) is one of the fibres. This projection is a Riemannian submersion 
onto V{M) equipped with the Wasserstein (see |Utj ). 

Consider the "horizontal" geodesies in D(M), which project to geodesies in V(M), 
see Fig.l. This projection means that instead of solutions of the Burgers equation 
()2.3|) we are interested only in the divergence of the corresponding fields, i.e. in 
how they act on densities. The latter is described by the infinitesimal version of the 
Monge-Ampere equation (j2.2j) . The flow of this infinitesimal version delivers a solu- 
tion to the Monge-Ampere equation for the potential of a gradient diffeomorphism, 
see |Brj . 

On the other hand, asymptotic directions on ST>(M) correspond to "almost ver- 
tical" geodesies in the space T>(M) with respect to the same flat L 2 metric. More 
precisely, an asymptotic geodesic (given by the Burgers equation) is "vertical" since 
it joins diffeomorphisms belonging to the same fiber ST>(M), i.e. it joins different 
pre-images of the point O on the quotient. Asymptotic directions are vertical initial 
vectors which correspond to geodesies having second order of tangency with the 
fiber. In other words, the projections of these geodesies to the base leave the initial 
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Figure 1. 71 is a geodesic in the space of densities V(M), 72 is 
a horizontal lift of 71 to P(M), 73 is an asymptotic (or 'vertical') 
geodesic in SD(M), 74 is an 'almost vertical' geodesic for SD(M) C 
T>(M) with the asymptotic direction X. 

point O very slowly. Indeed, recalling that taking the projection is the same as 
taking the gradient part, we identify here both of equations (jl.2|) : along with the 
Burgers equation they imply that divX = divX t = 0. These relations mean that 
the projection of the corresponding geodesic to the base V(M) is a curve which 
starts at the point O with zero velocity and zero acceleration, respectively. (For 
non-asymptotic directions only the velocity is zero, while the acceleration is not, 
since typical geodesies have the first order of tangency to SV{M).) 
The above discussion can be summarized in the following statement: 

Theorem 12. The Monge- Ampere equation in 2D fluid dynamics and (the 

infinitesimal form of) the Monge- Ampere equation in mass transport \2.2}) are pro- 
jections of the same dispersionless Burgers equation. 

Remark 13. Yet another appearance of the Monge- Ampere equation is related to the 
discrete version of the 2D Euler equation for a special energy functional, due to Moser 
and Veselov |MVj . A solution of the discretized Euler equation for an ideal fluid 
is a recursive sequence of diffeomorphisms. The Monge-Ampere equation provides 
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the constraint on the initial condition ensuring that all diffeomorphisms of the dis- 
crete geodesies are area-preserving. This restriction is similar to the Monge- Ampere 
condition singling out asymptotic directions among 2D divergence-free vector fields. 

2.3. Relations of diffeomorphism groups. One of the most interesting problems 
in the context of diffeomorphism groups is the existence of a shortest path or a 
geodesic between any two diffeomorphisms. This question is interesting already in 
dimension n = 1, for the group D(S l ) equipped with either L 2 or H 1 metric. Locally 
such a geodesic always exists (see |M2j ). 

In finite-dimensional geometry such questions are answered using the Hopf-Rinow 
theorem, but in infinite dimensions things become complicated. Grossman and Atkin 
constructed examples which show that even complete infinite-dimensional Hilbert 
Riemannian manifolds have points that cannot be joined by any geodesic. For 
instance, for the groups ST>(M n ) with n > 3 Shnirelman |Sh| constructed examples 
of pairs of diffeomorphisms for which a shortest path does not exist. 

In the study of this problem the following relation of the groups of diffeomorphisms 
of one- and two-dimensional manifolds can be useful. Consider the fibration of all 
diffeomorphisms of a surface M over the space of densities on M, which is exactly 
the context of optimal mass transport: 

SV{M) ^ V(M) -> C°°(M) . 

Here the left arrow is an isometric embedding with respect to the corresponding 
L 2 -metrics discussed above, and the right arrow is a Riemannian submersion into 
the space of densities C°°(M) equipped with the Wasserstein L 2 -metric. 

Now, consider a surface M 2 whose boundary is the circle S 1 . Similarly to the 
above, we have the following fibration 

SV (M 2 ) SV(M 2 ) -> V^S 1 ) 

where ST> (M 2 ) denotes the subgroup of smooth area-preserving diffeomorphisms 
of the surface that are pointwise fixed on the boundary S 1 . 

It would be interesting to study the Riemannian properties of this fibration. In 
this context the shortest path problem in the group ^(S 1 ) looks like the problem 
for an optimal path in the group ST>(M 2 ) or V(M 2 ), assuming that the projection 
to ©(S" 1 ) is a Riemannian submersion. The latter is a problem of optimal transport 
with a prescribed boundary map: we are connecting in an optimal way two disk 
diffeomorphisms which are liftings of the given diffeomorphisms of the circle. Here 
one could hope to employ the following reasoning. While the optimal map exists 
and is unique between any two convex domains, it automatically determines how 
the boundary is mapped. Hence one can almost never solve the same problem of 
finding the optimal map and simultaneously satisfy a particular boundary condition. 
This would allow one to conclude that the shortest path does not exist for almost 
all pairs of circle diffeomorphisms. 
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